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The microscopic formulas for the shear viscosity r/, the bulk viscosity and the 
corresponding relaxation times r-n- and rn of causal dissipative relativistic fluid- 
dynamics are obtained at finite temperature and chemical potential by using the 
Q . projection operator method. The non-triviality of the finite chemical potential cal- 
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culation is attributed to the arbitrariness of the operator definition for the bulk 
viscous pressure. We show that, when the operator definition for the bulk viscous 



> 

QQ ■ pressure 11 is appropriately chosen, the leading-order result of the ratio, Q over rn, 

^SJ , coincides with the same ratio obtained at vanishing chemical potential. We further 



discuss the physical meaning of the time-convolutionless (TCL) approximation to 
the memory function, which is adopted to derive the main formulas. We show that 
the TCL approximation violates the time reversal symmetry appropriately and leads 
results consistent with the quantum master equation obtained by van Hove. Further- 
more, this approximation can reproduce an exact relation for transport coefficients 
obtained by using the f-sum rule derived by Kadanoff and Martin. Our approach 
can reproduce also the result in Baier et al.(2008) Ref. ^ by taking into account 
the next-order correction to the TCL approximation, although this correction causes 
several problems. 



I. INTRODUCTION 



Relativistic hydrodynamics has been widely used to study various phenomena in many 
areas, such as heavy-ion collisions, relativistic astrophysics and cosmology. However, the 
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formulation of relativistic hydrodynamics is not trivial because the naive introduction of 
dissipations violates the relativistic causality, and leads to inconsistent behaviors in the 
property of stability jl, 2|. Thus the relativistic Navier-Stokes theory contains superluminal 
modes and is inadequate as a candidate of relativistic hydrodynamics 3|. 

Many different approaches have been studied for constructing the consistent relativistic 
hydrodynamics, including phenomenological approach j4-9|, kinetic approach 10, 11| and 



others 



12| . The common feature of these approaches is that the shear stress tensor and the 



bulk viscous pressure are treated as dynamical variables in addition to usual hydrodynamic 
variables (for example, energy density and fluid velocity), and the evolution of the shear 
stress tensor (the bulk viscous pressure) follows a dissipative equation characterized by 
the shear viscosity (the bulk viscosity) and the corresponding relaxation time. Different 
approach predicts different values for the relaxation times, and also gives different non- 
linear terms which can appear in the evolution equations as well. In the following, we call 
these theories the causal dissipative relativistic fluid-dynamics (CDRF) 

In hydrodynamics, all transport coefficients, such as the shear viscosity and the bulk vis- 
cosity, are inputs and should be determined from a microscopic theory or experiments. For 
the Navier-Stokes theory, the transport coefficients are usually calculated in two different ap- 
proaches. One is the kinetic approach based on the Boltzmann equation, and the other is the 
microscopic (field-theoretical) approach through the Green-Kubo-Nakano (GKN) formula. 
As known, the applicability of the Boltzmann equation is quite limited and the GKN formula 
are more general in application. However, it should be emphasized that both approaches 
give the same results in the classical dilute gas limit [isl [l^ . 

Differently from the Navier-Stokes theory, the computing methods for transport coeffi- 
cients in CDRF are not well established. In this work, we will apply the projection operator 
approach to derive the microscopic formulas for the transport coefficients in CDRF |l5l-ll8|. 
which can be viewed as generalization of the GKN formulas for the Navier-Stokes (Newto- 
nian) fluids. Although different approaches have been proposed [3,1^, 19|, as far as we know. 



this is the unique approach for which the consistency is confirmed by comparing to the 



results obtained from the Boltzmann equation 



18j. Moreover, it is verified that our formula 



satisfies the exact result obtained from the f-sum rule when applied to the (generalized) 



diffusion process 



17 



20 



2l|. 



The results at the vanishing chemical potential were already reported in Refs.|l6l. [ij 
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The purpose of this paper is to extend this analysis to the case of finite chemical potential. 

The non-triviality of the finite chemical potential calculation is attributed to the arbi- 
trariness of the operator definition of the bulk viscous pressure. At the vanishing chemical 
potential, this arbitrariness was removed by comparing with the result from the Boltzmann 
equation, but there is no corresponding result at finite chemical potential so far. We show 
that, if the bulk viscous pressure is defined appropriately, the leading-order result of the 
ratio of the bulk viscosity ( to the corresponding relaxation time rn coincides with the same 
ratio obtained at the vanishing chemical potential. 

Another purpose of this paper is to show the physical and mathematical meanings of 
the approximation used in the derivation. In the derivation of the transport coefficients of 
dissipative equations, we need to introduce non-trivial approximations to violate the time 
reversal symmetry. As far as we know, there is still no established method for this. In 
this paper, we use the so-called time-convolutionless (TCL) approximation. The meaning 
of this approximation was already discussed in Ref. In this paper, we argue that 

the quantum master equation obtained by van Hove and the f-sum rule satisfied in exact 
quantum time evolution processes are reproduced only when the TCL approximation to the 
memory function is applied. More importantly, in order to keep the consistency between 
the results of the Boltzmann equation and quantum field theory we have to use the TCL 
approximation, as was shown in Ref. 18|. A naive derivative expansion scheme does not 
satisfy these conditions. 

This paper is organized as follows. In Sec. [Tll the projection operator method is briefly 
summarized. In Sec. IIIIl the meaning of the TCL approximation is discussed. By using 
this approximation, the microscopic expressions for the shear viscosity, the bulk viscosity 
and the corresponding relaxation times are derived in Sec. IIVI Sec. |V] is devoted to the 
concluding remarks. We use natural units c = h = ks = 1- 



II. PROJECTION OPERATOR METHOD 
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-|2|. The 



Here we briefly summarize the projection operator method for CDRF 
difference between the calculations at vanishing and finite chemical potentials comes from 
the different definitions for the projection operators. 



Regardless of the presence of the chemical potential, the time evolution of an arbitrary 
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Heisenberg operator O is determined by the Heisenberg equation of motion, 

^0{t)=i[H,0{t)]^iLO{t), (2.1) 

where 0{t) = Qi^HtQ^-iHt ^^\^\^ }j ^j^g Hamiltonian and L is the Liouville operator. 

Let us introduce a set of operators for gross (in our case, the hydrodynamic) variables by 
A = {Aj}, z = 1, ■ ■ ■ ,77,. In order to implement coarse graining of the Heisenberg equation 



of motion, the projection operator is introduced. Following previous works, 16l-ll8l. l22|. we 
adopt the Mori projection operator, 

n 

P0= 5^(0,At)(A,AVA,. (2.2) 

The inner product here is given by Kubo's canonical correlation, 

(X, Y) = j^^ ^Tr[pe,e^^Xe-^^F], (2.3) 

where peq = e~^^/Tr[e~^^] with K = H — fiN and being a conserved charge. Note that 
the statistical expectation is taken not with H but with K. It is easy to check that 

{iLX,Y) = -{X,iLY), (2.4) 

even for the case of finite chemical potential. 

By using this projection operator, we can re-express the Heisenberg equation of motion 
for the hydrodynamic variables as follows, 

^A(t) = zA Ait) - 1^ rfrH(r)A(t - r) + m, (2-5) 

where A and H are (n x n) matrices and C, is an n-vector of operators. Their elements are 
given by 

zA,, = ^(zLA„4)(A,At)-i, (2.6) 

k 

S.,(t) = -9it)J2i^LQe^''^HLA,Al)iA,A^)~^, (2.7) 

k 

Ut) = Qe*^«*2LA, (2.8) 

where Q = 1 — P. The first and second terms on the right-hand side represent the collective 
oscillation and dissipation after coarse graining of time scale, respectively. The memory 
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function Hjj(t) can be expressed in terms of the time correlation of the third term, through 
the fluctuation-dissipation theorem. If the projection operator is chosen appropriately so 
as to collect all the macroscopic degrees of freedom, the third term will oscillate very fast 
and can be interpreted as the noise term. In the following calculation, the noise term is 
neglected. 

III. VIOLATION OF TIME REVERSAL SYMMETRY, SUM RULE AND 

COARSE GRAINING 

In this section we discuss the approximation which will be used in the following sections 
to derive our formulas of transport coefficients. See also the argument in Sec. V of Ref. fl?! . 

We have to keep in mind that the evolution equations of viscous fluids violate the time 
reversal symmetry, while the Heisenberg equation of motion (12. ip is time reversal symmetric. 
Thus, to derive the fluid dynamics from an underlying microscopic theory, some non-trivial 
operations are required. Notice that this is different from the case of the derivation of 
the fluid dynamics from the Boltzmann equation, because the Boltzmann equation already 
violates the time reversal symmetry by the assumption of the molecular chaos. In short, 
there are two steps to derive fluid dynamics from a microscopic theory: 1) the reduction of 
the number of the dynamical variables and 2) the violation of time reversal symmetry. In 
the Boltzmann equation approach, the second step is already done as was mentioned above 
and the first step corresponds to, for example, the Chapman-Enskog or moment expansion. 
Then the one-particle distribution function is replaced by hydrodynamic variables such as 
the energy density and the fiuid velocity. In the projection operator method, the first step 
corresponds to the appropriate choices of the projection operators (12. 2p . and the second step 
corresponds to the coarse graining for the memory function, which we discuss in this section. 

A. Violation of Time Reversal Symmetry 

As is commonly believed, the origin of the violation of the time reversal symmetry is 
attributed to the existence of two different time scales: one is macroscopic and the other is 
microscopic. The time reversal symmetry is violated when the degrees of freedom associated 
with the latter scale are neglected compared to the former. However, it is not easy to 
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implement this coarse graining systematically and there are several proposals 25|]. Here, we 



discuss the procedure introduced by van Hove in the derivation of the master equation from 



the Schrodinger equation 26[. In his derivation, the time variable t is rescaled as r which is 
defined by r = g^t with g being the coupling constant for the microscopic interaction. Then 
van Hove collected all terms which survive in the asymptotic limit of t — oo with r fixed 
and succeeded in deriving the quantum master equation. This limit is called the van Hove 
limit. Mathematically speaking, the van Hove limit corresponds to the procedure of picking 
up only secular terms which dominate the dynamics at the macroscopic time scale (which is 
realized in small coupling limit in his case), van Hove's argument indicates that to obtain 
macroscopic dissipative equation correctly, we have to collect only appropriate terms such 
as the secular terms. 

In the projection operator method, to violate the time reversal symmetry, we have to 
introduce a coarse graining of time scale in the memory function Hjj(t). As is discussed 



22l |. the memory function is given by the time correlation of the noise term, which is the 
consequence of the fluctuation-dissipation theorem. The time scale of the noise is much 
shorter than that of {Ai} and hence the time scale of '^ij{t) can be negligibly small comparing 
to that of {Ai} ^. Then the time-convolution integral of Eq. (12.51) is approximated as 

^A(t) = A{t) - dTS{r)A{t) + ^(t). (3.1) 



This is the TCL approximation ^. See also the discussion in Ref. [17[. On the other hand, if 
we carry out the Taylor expansion of A(t — r) in terms of r up to the next-to-leading order 
^, we have 

^A(t) = zA A(t) - drH(r)(A(t) - TdtA{t)) + ^(t). (3.2) 



^ If it is not the case, there are two possibihties. One is that the definition of the projection operator 
is stiU incomplete and {Ai} does not span the completely set of the gross variables. Then we have to 
generalize the definition of the projection operator. The other possibility is that there is no coarse-grained 
macroscopic theory and we have to solve the microscopic dynamics exactly. Thus the time dependence 
of the memory function can be the criterion to see whether the definition of the projection operator is 



appropriate or not. See also the discussion in Ref. [27 1. 
^ We do not call it the Markov approximation, because the memory effect can exist even after this approxi- 
mation. For example, the Maxwell-Cattaneo-Vernotte equation which we will discuss later has a memory 
effect even after the TCL approximation. This is also true for the derivation of transport coefficients of 
CDRF in the next section. 



This is different from the argument developed in Ref. [l^ ■ 
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Note that we still assume that the dominant contribution of the integral comes from r ~ 
and the upper limit of the integration is replaced by oo. 

At first glance, it might be considered that the latter approach (Eq. (13. 2p ) is more 
reliable because the next-order correction in the time-derivative expansion is considered. 
However, this is not trivial from the view point of the appropriate violation of the time 
reversal symmetry. Here we discuss which approximation is consistent with van Hove's 



argument 



26| by applying the projection operator method to derive the quantum master 



equation. The detailed derivation is shown in Appendix [Al As is shown by Eq. (lAlOp . 
van Hove's result is reproduced only when the TCL approximation is adopted. That is, 
the TCL approximation corresponds to the procedure of collecting all the secular terms in 
deriving the quantum master equation. On the other hand, when we consider the next-order 
correction to the TCL approximation, we pick up even irrelevant contributions. This is one 
of the evidences supporting that the TCL approximation may work better than Eq. (13. 2p 
in describing macroscopic physics. 



B. Exact Relation for Transport Coefficients 



Next, we show that the TCL approximation gives a consistent result with an exact relation 
obtained from the microscopic dynamics, while the next-order correction gives rise to an 



inconsistency |l7l. l20l. l21|. 



To show this, let us consider a non-relativistic diffusion process. Before deriving the 
(generalized) diffusion equation from the microscopic dynamics, we will point out that there 
is an exact relation for the dynamics of a conserved density following Refs. IT], |20|, |2l| . We 
consider the complex Schrodinger fields, ip and ip'^ , whose dynamics conserves the (spatial 
integration of the) number density defined by n = ip'^ip., as is shown in Ref. j20|. From the 
Noether's theorem, we can obtain the corresponding current operator J. Then there is the 
following relation. 



— V5=^(x-X')(n(x'))e,. 

m 



([J(x,t),n(x',t)])eg : 

This leads to the f-sum rule, 

— / — a;ImC^(k,a;) = -/3 / rf3x(J(x, 0), J(0, 0)) = -(n(x = 0))eg, 
y vr J m 



(3.3) 



(3.4) 
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where ( )eq denotes the thermal expectation value. Here we introduced the retarded Green 
function through 

-z([n(x,t),n(x',t')])e,^(t-tO = /^/ ^C«(k,c.)e-(*-*')e^'^(— (3.5) 

By using this relation, the exact form of the Laplace- Fourier transform of the time evolution 
of the conserved number density (5n(x, t) = n(x, t) — (n(x))eg is given by 

^ = -C^'ik, 0) - ^k^/? j rf='x(J(x, 0), J(0, 0)) + 0(1//), (3.6) 

where F(k) is an arbitrary function related to the initial condition. This is the exact relation 
which is obtained from quantum mechanics. Note that the above relation is independent of 
the form of the interaction term in the Hamiltonian. 



As was already emphasized in Refs. 



13, 



20 



2l| . if we assume the diffusion equation as 



the coarse-grained dynamics for 6n{x.), we cannot reproduce the exact result (13. 6p . In fact, 
the Laplace-Fourier transform of the diffusion equation is given by 

Sn^^{k,z) . C^(k,0) C^(k,0)Dk^ , 

-^^ = '^z z^ + ^'-'^ 

Differently from the exact result (13. 6p . one observes the following features: 1) there is a 1/2;^ 
term which does not disappear in Eq. (13.61) and 2) the coefficient for 1/z^ cannot reproduce 
the exact result, although it is not shown here. Thus the diffusion equation cannot reproduce 
the exact behavior obtained from quantum mechanics. 

However, if we consider the Maxwell-Cattaneo-Vernotte-type generalized diffusion equa- 
tion, 

dt5n{x, t) + VJ(x, t) = 0, (3.8) 
5t J(x, t) = -— Vn(x, t) - — J(x, t), (3.9) 

we can completely reproduce Eq. (l3.6p up to 0{l/z)^ as, 

^n^^(k,z) _ . C^(k,0) , . Dk^C^(O,0) , 

F(k) " ' i ' ^^-^^^ 

if the transport coefficients D and td satisfy the following relation, 

-r> = -C^) J ' '^""^ = /d3x(^n(x,0),MO,0)) - ^ 
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If we require our approximation to the projection operator method to be consistent with 
the f-sum rule, the derived D and To must satisfy Eq. fl3.1ip . In fact, as shown in Ref. j2o| . 
if define the projection operator with two gross variables, 6n{x.) and J(x), the equations of 
motion read 

dtSn{'k,t) +ik- 3ik,t) = 0, 

dt3{k,t) = -ikQ^6n{x,t) - [ c/rS^(k, r)J(k, t - r), (3.12) 

Jo 



where 



/c/3x(J(x,0),J(0,0)) 



/c/3x(5n(x,0),(5n(O,0))' ^ ' ^ 

When we adopt the TCL approximation to the memory function S-^(k, r), we find the 
following relations, 

- = / rfrS^(0,r), (3.14) 
td Jo 

— = n"". (3.15) 

To 

They coincide exactly with the exact result (13. 6p . 

On the other hand, if we consider the higher order correction as is done in Eq. (13. 2p . the 
definition oi I/td is modified as is shown in Eq. (65) of Ref. [l3], and one cannot reproduce 
the exact relation for the ratio D/td anymore. This is the second reason of why we do not 
take into account the next-order correction. 

The above argument is developed for specific examples, the quantum master equation and 
the generalized diffusion equation. Thus, exactly speaking, it is not obvious whether the 
next-order correction still causes inconsistency in the transport coefficients of CDRF. There 
may exist examples where the time-reversal symmetry can be violated by the derivative 
expansion used to derive Eq. (13. 2p . 

However, in the case of CDRF, there is another evidence to support the TCL approx- 



imation. The transport coefficients of CDRF can be calculatec 



even from the Boltzmann 



equation with the 14 moment approximation 10|. In Ref. 18|, it was shown that those 
results are consistent (not the same) with those from the projection operator when the TCL 
approximation is employed at vanishing chemical potential. The fact that the consistent 
results are obtained from two different approaches is very surprising and strongly suggests 
to employ the TCL approximation. 
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The expression of the relaxation time including the next-order correlation is explicitly 
given by Eq. (65) in Ref. 17|], and this expression is exactly the same as the result obtained 
in Ref. jsj. That is, our projection operator approach can reproduce the result of Ref. jsj 
by changing the approximation to the memory function. However, as was discussed in this 
section, we cannot obtain reasonable results if we include this next-order correction to the 
TCL approximation. This indicates that the derivative expansion, which is used in Ref. 
js], may not violate the time reversal symmetry appropriately. In fact, exactly speaking, 
hydrodynamics derived in this way corresponds to the relativistic Burnett equation and the 
corresponding transport coefficients are not those of CDRF, although it was assumed that 
they are equal in some literatures. Note that the Burnett equation has an intrinsic problem 



called the Bobylev instability and is essentially different from CDRF 36|, l37l| . Furthermore, 
the weak coupling limit of the result of Ref. sl was discussed in Ref. 19 1, but it is not 



consistent with the results from the Boltzmann equation with the 14 moment approximation 
as is discussed in Ref. |l8 |. 



IV. TRANSPORT COEFFICIENTS AT FINITE CHEMICAL POTENTIAL 



In this section, we will apply Eq. (12. 5 p and the TCL approximation to derive microscopic 
formulas for the shear viscosity, bulk viscosity and the corresponding relaxation times of 
CDRF. 



A. Shear viscosity and corresponding relaxation time 



Following Ref. 16|, let us consider a fluid flowing in the x direction with flnite velocity- 
gradient in the y direction. Thus the bulk viscous pressure does not show up. Then T"^ and 
T^^ are chosen as the gross variables, and the projection operator is given by 

PO = {0,T^''{-ky)) rj.Oxny\ I (0,T^''(-fc^)) rj.yxnV\ (A^\ 

{TO^{ky),T^^{-ky)) ^ ' {Ty^{ky),Tv^{-ky)) ^ 

Here we have implemented the Fourier transformation in space. Substituting it into Eq. 
(12. 5p . we obtain 

dtT'^^k^t) = -ikyTy%ky,t), (4.2) 

dtTy^{ky,t) = -ikyRiyT^^'ik^t) - [ dTE''{ky,T)Ty^{ky,t-T). (4.3) 

Jo 
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One can easily confirm that this equation is still symmetric under the time reversal op- 
eration, t ^ —t, by using the exact expression for S'^(/c^,t) 22|, |2J]. In order to violate 
the time reversal symmetry, we employ the coarse graining of time scale by using the TCL 
approximation. 



9fT°"(F,t) = -^A;^,T^'"(A;^^), 
dtTy^{k\t) = -ikyRlyT'>%ky,t) 



where 



l\,h.y 



{Ty^{ky),Ty^{-ky)) 



(4.4) 
(4.5) 

(4.6) 



As a result, the second term on the right hand side does not contain the time convolution 



17 



20 



22 



integral any more. See Refs. fl5 
memory function S'^(/c^, r) is given by 

lim^H'^^(/t^s 

with 



23| for details. The Laplace transform of the 



X^L{ky,s) 



X^^{ky,s) = / dte 
Jo 



■St 



{Ty%ky,t),Ty%-ky)) 



(4.7) 



(4.8) 



{Ty'={ky),Ty''{-ky)) 

On the other hand, the phenomenological equation for the shear stress tensor near the 
rest frame is 

r^^T^^(F) + T^^(F) = -r]{iky)u%ky), (4.9) 

where r) and are the shear viscosity and the corresponding relaxation time, respectively. 
By comparing Eq. (14.51) with Eq. (14.91) . we finally obtain the microscopic expressions for rj 
and r^- as 

1 



V 



lim -i?^,X"^(F,s) 



I3{e + P) s,ky-^o 13 

lim [ dtd^xe-'\Ty%x,t),Ty%0)) 



(4.10a) 



/3 / rf^x(T'^^(x),T'J^(0)) 



^ = lim lx'^\ky,s) 

j3 s,kv^o (5 



lim I dtd'^xe-'\Ty%x,t),Ty''(0)) 
^ / (i=^x(T^"(x),TJ'"(0)) 



(4.10b) 
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These expressions are formally the same as the previous results obtained at vanishing 



chemical potential 16|. The difference comes from the fact that, at finite chemical potential, 
the time evolution is governed by H but the ensemble average is taken with K = H — fiN. 
Thus it should be noted that the numerators of Eq. f l4.10p can be expressed using the 
analytic continuation of the Matsubara function only when T^^ and T^^ commute with 



In general, transport coefficients are related to the imaginary part of certain retarded 
Green functions. Thus one cannot obtain finite results for rj and unless the effect of 
interaction is considered. On the other hand, the dimensionless tj-Tt^ ratio defined hj ri/T.,r{e+ 
P) is finite even in the non-interacting case, because this ratio is given by the real part of the 
Green functions. Here, for the sake of simplicity, we consider the leading-order estimation, 
that is, we use the free gas approximation to calculate this ratio. 

Let us first consider a charged scalar boson described by the following Lagrangian, 

£ = d^(t)^d^'(t)-m^(t)'^(t). (4.11) 



As is discussed in Refs. 17|, l29|] , the canonical energy- momentum tensor for the scalar field 
theory is not well-defined. Instead, we use the improved energy-momentum tensor which 
reads, 

- g^"'{dp(p^d''<i)-m^(t)'^(P) 

- -{d^d" - g'"'d^)(P^(j). (4.12) 
3 

Then a straightforward calculation leads to 

rj _ J d3x(T?'^(x),T^^(0)) _ P 



T^{e + P) f d^^iT^^ (x) , (0) ) £ + P ' 
where the energy density e and the pressure P are given by 



(4.13) 



'^'P-E ( i + ^ ) (4 14) 



(27r)3 ^ yef^iEp^^i) - 1 e/3(i?p+M) _ i 

P = I d'P If I + I ) (415) 

' (27r)3Ep 3 \ePiEp~f^) - 1 e/3(^p+M) - 1 ' ' ^ ' 



respectively. Here Ep = ^/p'^ + w?. In this derivation, the temperature-independent diver- 
gent term is neglected. 
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FIG. 1: The tj-Ttt ratio, 7j/Tn{£ + P), as a function of temperature and chemical potential. 

The temperature and chemical potential dependence of the //-r^ ratio given by Eq. f l4.13p 
is shown in Fig. [T] We use the pion mass m = 140 MeV. One can see that the ratio 
monotonically increases as a function of temperature, and converges to = 0.25 which is 
the result for the massless bosons. 

Compared to the temperature dependence, the chemical potential dependence is not easily 
recognized. In order to see the chemical potential dependence more clearly, we plotted the 
temperature dependence at fixed chemical potentials, /i = 50, 100 and 130 MeV in Fig. [21 
One can see that the weak chemical potential dependence still exists and the ?7-r,r ratio 
decreases slowly as the chemical potential grows. 



As is discussed in Refs. j2|, this ratio is directly related to the propagation speed of 
signals in CDRF. Thus, for the theory to be relativistically causal, the //-r^ ratio should not 
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FIG. 2: The r/-T^ ratio as a function of temperature at fixed chemical potentials. The solid, dashed 
and dot-dashed lines represent the ratio at /u = 50, 100 and 130 MeV, respectively. 

be larger than one. This condition is satisfied for any temperature and chemical potential 
for our calculation. 

These results will be modified by the interaction. As a matter of fact there are attempts 
to estimate this ratio including the effect of interaction by using the lattice QCD simulation, 
see, for example, Ref. [30 1. 



B. Bulk viscosity and corresponding relaxation time 



Following Ref. 



171], we consider a perturbation in an infinite fluid in thermal equilib- 



rium having a planar symmetry in the (?/, z) plane. All the quantities associated with the 
perturbed fluid dynamics vary spatially only along the x direction. In this case, the fluid 
velocity points to the x direction. Then and 11 are chosen as the gross variables and the 
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projection operator is given by 

PO 



(TO^(A;^),TO^(-A;^)) 



(4.16) 



(n(fc^),n(-fc^)) 

Here the operator of the bulk viscous pressure 11 is defined by the deviation from the equi- 
hbrium pressure. 

Differently from the case of the shear stress tensor, the definition of the operator of the 
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ail, 



bulk viscous pressure is changed by finite chemical potential. As is discussed in Ref. 
the transport coefficients are expressed by the vanishing momentum limit of commutators. 
Thus even if we added a term which commutes with the Hamiltonian of our system to 
the definition of 11, the final result is not affected except for the arbitrariness related to 
renormalization. In fact, from the view point of renormalization, we should add appropriate 



operators to the original definition of 11 3ll]. For the case of the vanishing chemical potential, 
we used 

1 ^ fip 

n = - ^(r^ - (r%,) - — (r°° - (T°°)e,). (4.17) 



i=l 



Note that the last term vanishes in the massless limit and does not violate the conformal 
property. The same modification of the bulk viscosity was proposed in Ref. 



31 



32|. As the 



justification of this modification, we would like to point out that the bulk viscosity calculated 
with our microscopic formula is consistent with the result of the Boltzmann equation only 
when we use Eq. fl4.17p as the operator definition of the bulk viscous pressure. 

At finite chemical potential, the operator of the conserved number density n commutes 
with the Hamiltonian, and the definition of H can be modified as 

1 



fdP 

n = -(T- - (T-)e,) - (^^ 



dP\ 



Meg)- 



(4.18) 



Note that (dP/dn)^ vanishes in the massless limit. See Appendix IB] for details. Thus this 
new definition does not violate the conformal proper ty of the bulk viscous pressure. The 
same modification of H was discussed in Refs. 
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32| 



Substituting Eq. fl4.16p into Eq. 02.51) . we obtain 



dtT"%k'',t) = -zA;"H(r,t), 



dtU{k\t) = -ik''R^.T°''{k\t) - [ rfrS^(A;^t)H(A;^^ 

Jo 



(4.19) 
(4.20) 
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where 

n _ (n(fc-),n(-fc-)) 



Similarly to the previous subsection, the TCL approximation is already employed [17|. The 
Laplace transformation of the memory function is given by 

lim = -^r^^^ -, (4.22) 



where 







On the other hand, the phenomenological equation of the bulk viscous pressure is given 

by 

rn^n(A:^) + n(F) = -C{i}z^)u\k^), (4.24) 

where C and rn are the bulk viscosity and the corresponding relaxation time, respectively. 

The microscopic expressions for and rn are obtained by comparing Eq. (14.201) with Eq. 
(I4.24p , employing the coarse graining of time jl3] , 



lim -R%X'''^{k\s) 



I3{e + P) s,k-^o (3 

lim ! c/trf3xe^'*(n(x,t),n(0)) 

/3 j c/3x(T°^(x),T°^(0)) 

— = lim -X^^ik\s) 



(4.25a) 



lim / dtd''xe-'*(n(x,t),n(0)) 

p [ rf^x(n(x),n(o)) 



(4.25b) 



Similarly to the case of the shear viscosity, we calculate the (-tu ratio defined by (/tyi{b + 
P) by applying our results to a non-interacting charged scalar boson. By using Eqs. f l4.12p 
and (I4.18p . we can find 11. Then the ^-Tn ratio is calculated as 

C ^ /^Mn(x),n(o)) 

Tn{e + P) /#x(TO-(x),ro-'(0)) 
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where the sound velocity Cg is defined by 



.^^m .(^ji) ,4.7) 



\de J s_ \de J ^ e + P \dn 

n 

where e and P are defined in Refs. fl4.14p and f l4.15p . respectively, and 

n=[^( I I ^ (4 28) 

In this derivation, we neglected the temperature independent divergent term. One can easily 
check that the ^-rn ratio disappears in the massless limit. 

Differently from the case of the vanishing chemical potential, this ratio at the finite 
chemical potential has not yet been calculated in the frame work of coupled Boltzmann 
equation with particles and anti-particles. Thus we cannot discuss the validity of this result 
from the view point of consistency with the kinetic theory. However, we find that the 
functional form of this ratio is completely equivalent to the result for the vanishing chemical 
potential shown in Ref . l3] • This is not observed if we ignore the last term of Eq. f l4.18p 
which is newly added. Furthermore, the functional form of the 77- ratio is not changed 
by the introduction of the chemical potential as was shown in Eq. ( 14.13^ . Because of these 
facts, we believe that our result is reasonable and the bulk viscous pressure should be defined 
by Eq. f l4.18p at finite chemical potential. 

The temperature and chemical potential dependence of the C-rn ratio given by Eq. f l4.26p 
is shown in Fig. |3l We used the pion mass m = 140 MeV. One can see that the C-rn ratio 
monotonically decreases as a function of temperature, and finally vanishes. This is because, 
at very high temperature, the existence of mass is negligible and then the system effectively 
restores the conformal symmetry. 

In Fig. m the temperature dependence of the C-tji ratio is plotted at fixed chemical 
potentials, /i = 0, 50 and 130 MeV. At /i = 0, the ratio exhibits maximum at the vanishing 
chemical potential. However, it decreases more quickly for smaller chemical potentials. Then 
at higher temperature, this ratio becomes larger as the chemical potential increases. 

Similarly to the case of the shear viscosity, this ratio is directly related to the propagation 
speed of signals in CDRF. Thus, for the theory to be relativistically causal, the ratio should 
not be larger than one. This condition is satisfied for any temperature and chemical potential 
for our results. 
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FIG. 3: The C-^n ratio, C/7"n(e + -P), as a function of temperature and chemical potential. 

C. fermion 



So far, we have discussed the non-interacting charged scalar bosons. Parallelly, we can 
apply the same calculations to non-interacting fermions. 
Then we obtain 

^ 0, (4.29a) 



C _ {l-cl){e + P)-l{e-3P) 



(4.29b) 



Tn{e + P) e + P 

There results are the same as those obtained for the vanishing chemical potential case 
. The reason for the vanishing of ti/tt^^s + P) for fermions is because that the contribu- 
tion to this quantity from the pair annihilation and creation (PAC) processes cancels the 
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FIG. 4: The C-''"n ratio as a function of temperature at fixed chemical potentials. The solid, dashed 
and dot-dashed lines represent the ratio at /U = 0, 50 and 130 MeV, respectively. 

contribution from non-PAC pro cesses. More discussion can be found in Ref. 

As is emphasized in Ref. 18|, these calculations are performed at leading order and will 



be modified by the effect of interactions. For example, the exact expressions for t^^/P and 
rn//3 are given by the ratio of the real and imaginary parts of the retarded Green's function 
of T^^ and 11, respectively 17j. To leading order, the real part is given by the result of 
the free-gas approximation, while the imaginary part is not. That is, the leading-order 
calculation violates this exact relation and may lead to inconsistent results. The vanishing 
ri/{TT^{e + P)) for fermions could also be rendered finite by a more complete calculation 
including interactions. 

In addition, usually fermions interact by exchanging bosons. In such a mixed system of 
bosons and fermions, ?7/(r^(£ + P)) is given by the contributions from both fermions and 
bosons and takes a finite value even if the contribution from fermions vanishes. 
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V. CONCLUDING REMARKS 

In this paper, the microscopic formulas for the shear viscosity rj, the bulk viscosity (, and 
corresponding relaxation times r,r and rn of causal dissipative relativistic fluid-dynamics are 
obtained at finite temperature and chemical potential. 

Before obtaining these formulas, we first discussed the theoretical and mathematical 
meanings of the TCL approximation, which is used to violate the time reversal symmetry 
possessed by microscopic dynamics. To examine the validity of the TCL approximation, we 
applied the projection operator method to the derivation of the quantum master equation 
and to the generalized diffusion equations. By using the TCL approximation, we could collect 
only the secular terms appropriately and succeeded in deriving the quantum master equation 
obtained by van Hove. This means that the TCL approximation can pick up appropriate 
secular terms which should be kept for deriving dissipative equations. When the TCL 
approximation is applied to derive a coarse-grained equation of the non-relativistic diffusion 
process, the diffusion coefficient and the corresponding relaxation time satisfy the exact 
relation, the f-sum rule, which is obtained from quantum mechanics. On the other hand, 
we cannot reproduce these results, if we consider the next-order time-derivative correction 
to the TCL approximation. 



Moreover, as is discussed in Ref. 18|, the shear and bulk viscosities and corresponding 
relaxation times obtained by using the projection operator method can be consistent with 
the results from the Boltzmann equation only when the TCL approximation is applied. Thus 
we conclude that, to violate the time reversal symmetry appropriately, we should apply the 
TCL approximation and should not consider the next-order correction. 

However it is worth mentioning that, if we take into account this next-order correction 



to the TCL approximation, we obtain the result given by Eq. (65) in Ref. 17||, which is the 



same formula obtained in Ref. 
shown in Ref. 
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The quantitative difference of our result and theirs are 
18[ for the vanishing chemical potential. 
Next, we derived the formulas of the transport coefficients of the causal dissipative rel- 
ativistic fluid dynamics at finite temperature and chemical potential by using the TCL 
approximated projection operator method. The formulas for the shear and bulk viscosities 
are given by Eqs. f l4.10p and fl4.25p . respectively. 

In the calculation of the bulk viscosity, the operator expression of the bulk viscous pressure 
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should be chosen appropriately. In this work, we applied Eq. fl4.18p . the last term of which 
does not exist in the calculation at the vanishing chemical potential fl?! . 

The transport coefficients r], (, t^^ and rn are not finite in the free-gas approximation 
because they are proportional to the imaginary part of the Green functions. However, the 
ratios ///r^ and C/'^n ^''^^ still finite because those are calculated from the real parts of the 
Green functions. Then we obtain 

— = (3-a)P, (5.1) 



where a = 2 for the charged scalar boson and a = 3 for fermion. The functional forms 
of these ratios are completely the same as the results for the vanishing chemical potential 
shown in Ref . [ig], [i3| . This equivalence is not observed if we ignore the last term of Eq. 
( ]4.18p which is newly added. It should be noted that these expressions are confirmed only 
for the leading order calculations and it has not yet known how these relations are modified 
under the effects of interaction. 

This may suggest the consistency of our calculations. However, to confirm it more pre- 
cisely, we should compare our results with those from the Boltzmann equation including 
particles and anti-particles. In the case of the vanishing chemical potential, the same ratios 
were calculated from the simple Boltzmann equation and confirmed that these are consistent 
when quantum correlations which are not included in the Boltzmann equation are neglected 

n 

[18l |. The consistency check for the case of finite chemical potential is left for future task. 

Besides the viscosities and corresponding relaxation times, CDRF can contain other trans- 
port coefficients. For example, the heat conductivity and the corresponding relaxation time 
are other important transport coefficients ^. These are under investigation. 
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^ As a matter of fact, the moment expansion of the Boltzmann equation predicts a lot of non-linear terms. 
However, we have to carefully pick up only terms which arc consistent from the viewpoint of the order of 
the Knudsen number 33 1. 
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Appendix A: The van Hove limit and the TCL approximation 



In order to see the relation between the van Hove hmit and the TCL approximation, we 



reformulate the result of van 



3ove 



by using the projection operator method, following 



the argument of Nakajima 34] and Zwanzig 351 ] . 

We consider a quantum system whose Hamiltonian is given by 

H = Ho + gV, 



(Al) 



where Hq is the non-perturbed Hamiltonian and the Fock space is defined by the eigenstates 
of Ho, 



Ho\Ea) = E^\Ea), {Ea\E'a') = 6{Ea - E'^,)5{a - a'), 



(A2) 



where a is a quantum number except for energy. Without loss of generalit y, w e assume that 
the interaction gV contains only off-diagonal parts in this representation 26]. 
The dynamics of the density matrix is given by 



dtp{t) = -iLp{t). 



(A3) 



This equation is re-expressed as the following form by introducing the general projection 
operator P satisfying = P, 



dtPpit) = -iPLPpit) - iPLe-^'^'Qp{0) + iPL [ dse-^'^'QiLP p{t - s) 

Jo 



(A4) 



where Q = 1 — P. This is easily derived by using the same method used in deriving Eq. 
(12 .Sp or using the following operator identity. 



-iLt 



Pe-'^' + e-''^^'Q - I [ dse-^'^'QLPe 

Jo 



■iL{t-s) 



(A5) 



To reproduce the master equation obtained by van Hove, following Ref. [3J], we specify 
P as a projection operator which extracts the diagonal parts of an arbitrary operator A, 



a\PA\E'^„ a') = {E^, a\A\E'^,, a')5{E^ - E'^,)5{a - a'). 



(A6) 



Note that the projection operator defined above is a superoperator. Following Ref. 26], 
we choose the initial condition so that p(0) does not contain off-diagonal components, i.e., 
Qp{0) = 0. Then Eq. becomes 



dtpnit) = -iPLpoit) + iPL [ dse-'^^'^'QiLpDit - s), 

Jo 



(AT) 
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where poit) = Pp(t) is the diagonal part of p(t). 

Corresponding to the Hamiltonian flAll) . the operator L is separated as 

L = Lo + Li, (A8) 

where Lq and Lj are the Liouville operators with Hq and gV, respectively. One can easily 
check that PLPA = and PLqA = because of the definition of P. Then, from Eq. (lA7p . 
the equation for poit) is re-expressed as 

dtpoit) = PiLi [ dse~^'^'QiLipD{t -s)^ PiLj [ dse-^'^'^'QiLipoit - s). (A9) 
Jo Jo 

Note that, to compare the result of van Hove, it is enough to expand e~*^*^^ in terms of the 
interaction strength g and keep only the lowest order term (more exactly, these higher order 
terms disappear in the van Hove limit). About this expansion, see, for example, Ref. 23 1. 

Let us introduce a probability ^^(a, t) = {Ea\p{t)\Ea) for a particle to be found at state 
a at time t. Then the equation for this probability is 

dtPE{a,t) = - j^ds J dE'^,da'2g^\{Ea\V\E'a')\^cos{E'^,-E^)s{PEia,t-s)-PE'ia',t- 
ds J dE'^,da'2g^\{Ea\V\E'a')\'^ cos{E'^, - E^)s {PEia,t) - PE'ia',t)} 
dE'^,da'2g'\{Ea\V\E'a')\' ( lim MElJz3iK] {Pj,{a,t) - Pe^^' ,t)} 

U'^oo E', - Ea J 



dE'^,da'w{a] a') {PE{a', t) — Pe{c(, t)} , (i 



where 



w{a]a') = 27ig^\{Ea\V\Ea')\'^6{Ea, - E'^,). (All) 

From the first to the second line on the r.h.s., we used the TCL approximation by replacing 
the upper limit of the integration from t to oo, and the index of Pe from t — s to t. The 
last result is completely equivalent to the master equation obtained by van Hove by using 
a different method j26|. As is well known, the transition matrix w coincides with Fermi's 
golden rule. 

As was mentioned in the text, van Hove obtained this result by collecting all the secular 
terms under the van Hove limit. In the projection operator method, the same result is 
reproduced using the TCL approximation. That is, the TCL approximation plays the same 
role as the van Hove limit and collects only relevant terms in the coarse graining of time 
scale. However, if we consider the next-order correction to the TCL approximation in Eq. 



(lAlOp . we cannot reproduce van Hove's result. 
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Appendix B: thermodynamic relations 



We will show that {dP/dn)^ vanishes in the massless limit. 
Note that 

= sdT + ndii. 



(Bl) 



From this relation, we obtain 

1 



{dP\ 
\dn)e 



Det 



dn 

d^l y rp 



T 









in 



where 



Det = if) 



dn\ ^ ( ds 



dn 
dT 



f dn 



T (— 



Thus {dP/dn)s is expressed as 



{dn\ 



T 



'de\ fde 



dT 
1 



Det 



(B2) 



(B3) 



(B4) 



In the conformal limit where e = 3P, {de/dfi)rp = 3n and (de/dT)^ = 3s. Thus (dP/dn) 
vanishes. 
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